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ABSTRACT 



Screw instability of the magnetic field connecting a rotating black hole (BH) with its 
surrounding disk is discussed based on the model of the coexistence of the Blandford-Znajek 
(BZ) process and the magnetic coupling (MC) process (CEBZMC). A criterion for the screw 
instability with the state of CEBZMC is derived based on the calculations of the poloidal 
1/^ \ and toroidal components of the magnetic field on the disk. It is shown by the criterion that 

^ \ the screw instability will occur, if the BH spin and the power-law index for the variation 
Q ■ of the magnetic field on the disk are greater than some critical values. It turns out that 
. the instability occurs outside some critical radii on the disk. It is argued that the state of 
^ \ CEBZMC always accompanies the screw instability. In addtition, we show that the screw 

' instability contributes only a small fraction of magnetic extraction of energy from a rotating 

r^: BH. 
o 

c/3 . Subject headings: accretion, accretion disk — black hole physics — magnetic fields — insta- 

bility 



■ 1. INTRODUCTION 

It is well known that the magnetic field configurations with both poloidal and toroidal components 
can be screw instable (Kadomtsev 1966; Bateman 1978). According to the Kruskal-Shafranov criterion 
(Kadomtsev 1966), the screw instability will occur, if the toroidal magnetic field becomes so strong 
that the magnetic field line turns around itself about once. Recently some authors discussed the screw 
instability in black hole (BH) magnetosphere. Gruzinov (1999) argued that a Kerr BH, being connected 
with a disk by a bunch of closed field hues, can flare quasi-periodically, and this process can be regarded 
as a new mechanism for extracting rotational energy from the BH. Li (2000a) discussed the screw 
instability of the magnetic field in the Blandford-Znajek (BZ) process (Blandford & Znajek 1977), which 
results in a stringent upper bound to the BZ power. Tomimatsu et el. (2001) studied the condition of 
the screw instability in the framework of the variation principle, and argued that the field-line rotation 
has a stabilizing effect against the screw instability. 
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Recently much attention has been paid on the magnetic couphng (MC) process, where energy and 
angular momentum are extracted from a fast-rotating BH to its surrounding disk by virtue of closed 
magnetic field lines (Blandford 1999; Li 2000b, 2002a; Wang, Xiao & Lei 2002, hereafter WXL; Wang, 
Lei & Ma 2003, hereafter WLM). The MC process can be used to explain a very steep emissivity in the 
inner region of the disk, which is found by the recent XMM-Newton observation of the nearby bright 
Seyfert 1 galaxy MCG-6-30-15 (Wilms et al. 2001; Li 2002b; WLM). Compared with the BZ process, the 
disk load in the MC process is much better understood than the remote load. Very recently we discussed 
the condition for the coexistence of the BZ and MC processes (CEBZMC), and some electromagnetic 
quantities, such as electric field, magnetic field and electric current in the BH magnetosphere, are 
determined (Wang et al. 2003, hereafter WMLY). 

The above works motivate us to discuss the screw instability in CEBZMC based on the Kruskal- 
Shafranov criterion and our calculations of the poloidal and toroidal components of the magnetic field. 
This paper is organized as follows. In §2 we give a brief description on our model of CEBZMC. In §3 the 
criterion of the screw instability is derived based on the Kruskal-Shafranov criterion and our calculations 
of the poloidal and toroidal components of the magnetic field on the disk. It turns out that the screw 
instability is related intimately to the two parameters, i.e., the BH spin and the power-law index for 
the variation of the magnetic field on the disk. In addition, the correlation of the screw instability with 
CEBZMC is discussed. It is shown that the state of CEBZMC always accompanies the screw instability, 
and the region of the screw instability can be determined by the criterion for the given values of the BH 
spin and the power-law index. In §4 we discuss the effect of the screw instability on magnetic extraction 
of energy from the rotating BH. It is shown that the screw instability contributes only a small fraction 
of magnetic extraction of energy from a rotating BH. Finally, in §5, we summarize our main results. 

Throughout this paper the geometric units G — c— 1 are used. 



2. DESCRIPTION OF THE MODEL OF CEBZMC 

Considering that the remote load in the BZ process and the load disk in the MC process are 
connected with a rotating BH by open and closed magnetic field lines, respectively, we studied a model 
of CEBZMC (WXL, WMLY). The poloidal configuration of the magnetic field is shown in Figure 1, 
where rj„ and Tout are the radii of the inner and outer boundaries of the MC region, respectively. The 
angle 6m indicates the angular boundary between the open and closed field lines on the horizon, and 
is the lower boundary angle for the closed field lines. The angles 9m and 9l are connected with rout 
and rin by the highest and the lowest closed field lines, respectively. Throughout this paper 6l = 0.457r 
is assumed in calculation. Following Blandford (1976), we assume that the poloidal magnetic field 
on the disk varies with the radial coordinate as a power law, 

B'n oc r", (1) 

Hereafter the subscript is used to indicate the quantities on the disk. The parameter n is the 
power-law index and ^ = r/rms is the radial coordinate on the disk, which is defined in terms of the 
radius Vms = MXms '^^ ^^e marginally stable orbit (Novikov & Thorne 1973). 

An equivalent circuit for the BZ and MC processes is used as shown in Figure 2, in which the 
symbols have the same meanings as given in WXL. 
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A stationary, axisymmetric magnetosphere anchored in a Kerr BH and its surrounding disk can 
be described in Boyer-Lindquist coordinates, and the concerned metric parameters are given as follows 
(MacDonald and Thorne 1982, hereafter MT). 



S2 = (r2 + a^f - a^A sin^ 9, = + a" cos^ 
A = + - 2Mr, w = (S/p) sin 9. 



(2) 



where M and a = J/M are the mass and the specific angular momentum of a Kerr BH, respectively. 
The current in each loop of the equivalent circuit for the BZ and MC processes is expressed by 

' + (^•,,/2.)^. (3) 



where the quantities in equation (3) have the same meaning as given in WXL. In WMLY we derived 
the poloidal currents on the horizon as follows. 

^H,BzMn)^Io^^^^^^y 0<9<9m, (4) 

^kMci^*,o,n) = io ^Jjl^J^^l^y eM<e<9,, (5) 

lo = B^^M ^ B4 (M/Mq) X 1.48 X 10^° A, (6) 

where the subscripts "BZ" and "MC" represent the quantities involving the BZ and MC processes, 
respectively, and B4 is the strength of the magnetic field in units of 10^ G. In equations (4) and (5) the 
parameter a* = a/M is the BH spin, and the parameter q is defined as q — y^l — a^. The quantity Bjj 
is the poloidal component of the magnetic field on the horizon. Hereafter the subscript "if" is used to 
indicate the quantities on the horizon. 

The parameters k and f3 are the ratios of the angular velocity of the open and closed magnetic 
field lines to that of the BH, respectively. Usually, k = 0.5 is taken for the optimal BZ power, while (3 
depends on the BH spin a* and ^. Since the closed field lines are frozen in the disk, P is expressed by 



+ a* 



(7) 



where fin and flo are the angular velocities of the BH and the disk, respectively. 



The poloidal current flowing on the disk is equal to that flowing on the horizon in the same loop 
of Figure 2, i.e.. 
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where ^ is related to the angular coordinate 6 on the horizon by the mapping relation given in WMLY, 



cos9 — cos9l — G {a*;^,n) d^, (9) 



where 



2 V(l + ^Ix^ni + 2a2x;^6) (1 _ 2^-2^-1 + a^x-te-^) 

The condition for CEBZMC is discussed in WMLY based on the following assumptions: 

(1) The theory of a stationary, axisymmetric magnetosphere anchored in the BH and its surrounding 
disk formulated in MT is applicable not only to the BZ process but also to the MC process. The 
magnetosphere is assumed to be force-free outside the BH and the disk. 

(2) The disk is both stable and perfectly conducting, and the closed magnetic field lines are frozen 
in the disk. The disk is thin and Keplerian, lies in the equatorial plane of the BH with the inner 
boundary being at the marginally stable orbit. 

(3) The magnetic field is assumed to be constant on the horizon, and to vary as a power law with 
the radial coordinate of the disk as given by equation (1). 

(4) The magnetic flux connecting a BH with its surrounding disk takes precedence over that con- 
necting the BH with the remote load. 

It is found that the state of CEBZMC will occTir, provided that the BH spin a* and the power-law 
index n are greater than some critical values. In the following sections we adopt the above assumptions 
to discuss the screw instability. 



3. CRITERIOR OF SCREW INSTABILITY IN CEBZMC 

3.1. Derivation of Criterion 

According to the the Kruskal-Shafranov criterion, the screw instability will occur, if the magnetic 
field line turns around itself about once, i.e., 

{271WD/L) B^j,/Bl < 1, (11) 

where and Br^ are the poloidal and toroidal components of the magnetic field on the disk, respec- 
tively, and ro£) is the cylindrical radius on the disk and reads 

= T^nlpD = eMxLyn^^!Px^i+2^!Px^- (12) 
In deriving equation (12) the Kerr metric parameters given in equation (2) are used. The quantity L 
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in equation (11) is the poloidal length of the closed field line connecting the BH with the disk. The 
criterion (11) can be rewritten as 

Bl/Bl = {Bl/B^^) {B'^Bl) < L/{27rzun). (13) 

Considering the continuum of magnetic flux between the two adjacent magnetic surfaces connecting 
the BH and the disk, we have 

S^27r {vjp\^^^ d9 = -S^27r (vjp / Va) dr. (14) 

Incorporating equation (14) with the mapping relation in WMLY, we express the ratio B^/B^ as 

B'n 2(l + g)v/l + a2e-2x-t-2e-ix.-' 



-G(a.;e,n). (15) 



The ratio B^/Bjj can be calculated by using equation (8), and the radial current density flowing 
on the disk can be written as 

fo (a*; n) = I^^MC ^)/ (STTOTij). (16) 

By using Ampere's law the toroidal magnetic field on the disk is expressed by 

= 47rj^ X n, (17) 
where n is the unit vector normal to the disk. Incorporating equations (8), (12), (16) and (17), we have 



B'h _ exL yi + aU-'x;ni + 2al^-'xJs [2 csc^ ^ - (1 - q)] 
Bl 2a, (1 - 13) 

Substituting equations (15) and (18) into the criterion (13), we have 



(18) 



{2T:wD/L)F{a,-i,n) <l, (19) 

where 



Since the poloidal length of the closed field line L can be estimated approximately as half of the 
circumference with diameter D, we have 
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L = AD = A 



r {grrf'dr, (21) 

J Th 



where D is the proper distance from the horizon to the place of disk where the closed field line penetrates, 
and A O.Stt is taken. The metric parameter g„ is expressed as 

= P = r- + a--2Mr^ ^xL + af- 2CxL ' ^^^^ 
Incorporating equations (12), (21) and (22), we have 

L/{27rwn) = ^/^"^ T (23) 



where = ru/rms = (1 + l)/x^s- Therefore equations (19), (20) and (23) consist of the criterion of 
the screw instability in our model. 

Inspecting the criterion (11), we find that the screw instability may occur, if the ratio of 5^ to 
is small enough. We can give an analysis before calculation. In fact, we find that Bj) is proportional 
to approximately by combining equations (6), (8), (16) and (17). Thus the ratio B^^f B^ should 
be proportional to by considering equation (1), and it is exactly the leading factor of the function 
F {a^; ^,n). In addition, from equations (8), (16) and (17) we find that the larger value of the BH spin 
results in the larger values of and B]j. Therefore we have the following conjectures on the screw 
instability. 

(1) The screw instability would occur at some place far away from the inner edge of the disk. 

(2) The screw instability would occur, provided that the BH spin a* and the power-law index n are 
greater than some critical values. 

The above conjectures are verified by the calculations on the criterion (19) as shown in Figure 3. 

In Figure 3a {2tiwbI L) F (a*; n) < 1 holds for 0.8903 < a* < 1 with n = 4, and in Figure 3b it 
holds for n > 3.34 with a* = 0.998. These results imply that the screw instability occurs, if the BH 
spin and the power-law index are greater than some critical values. 



3.2. Screw Instability and CEBZMC 

In WMLY we argued that the state of CEBZMC will occur, if the parameters a* and n are greater 
than some critical values. So it is tempting for us to discuss the correlation of the screw instabihty 
with CEBZMC. Based on WMLY and the criterion (19) we have the critical fine for CEBZMC and the 
critical lines for the screw instability in the parameter space consisting of a* and n as shown in Figure 
4. 

From Figure 4 we obtain the following results. 
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(1) There is very little difference among the critical lines for the screw instability with different 
values of the factor A, and these critical lines can not be distinguished in Figure 4a, and are very near to 
the critical line for CEBZMC (thick solid line) as shown in Figure 4b. Therefore, the screw instability 
is rather insensitive to the values of the factor A. 

(2) It is shown in Figure 4b that the critical line for CEBZMC corresponds to the larger values of 
a* and n than the critical lines for the screw instability. It implies that the state of CEBZMC always 
accompanies the screw instability. 



3.3. Disk Region for Screw Instability 

From the above discussion we conclude that the minimum radial coordinate ^screw for the screw 
instability can be determined by the following equation, 

{^-KWd/L) F (a*; Cscre«;, ^) = 1- (24) 

Equation (24) imphes that the disk region for the screw instability is 

^screw <^< 00. (25) 

By using the mapping relation (9) and equation (25), we obtain the corresponding angular region 
on the horizon as follows, 

Om < < screw (26) 

By using equations (19), (20) and (23) we have the contours of constant values of ^screw for the 
screw instability in the parameter space consisting of and n as shown in Figure 5. 

From Figure 5 we find that the screw instability could even take place in the inner region of the 
disk with a small value of iscrew) provided that a* and n are great enough. Since the large values of a* 
and n in the MC process are supported by the recent XMM-Newton observation for producing a very 
steep emissivity in the inner region of the disk (WLM, WMLY), these values imply not only CEBZMC, 
but also the screw instability. 



4. SCREW INSTABILITY AND ENERGY EXTRACTION FROM A ROTATING BH 

Gruzinov (1999) argued that the screw instability can be regarded as a new mechanism of extracting 
rotational energy from a BH. His argument is based on a model with a bunch of closed magnetic field 
lines connecting a Kerr BH with a disk. Compared with Gruzinov's model, the magnetic field in our 
model is axisymmetric, and the region of the screw instability is determined by equation (25) rather 
than is given at random. 

Considering the energy released due to the screw instability, the total power of the magnetic 
extraction from the BH consists of the following three parts, i.e.. 
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total 



Pbz + Pmc + P screw ■> 

(27) 



where Pbz-, Pmc and Pscrew are the BZ power, the MC power and the average power due to the screw 
instabihty, respectively. Incorporating the region for screw instabihty given by equation (26) with the 
expressions for the BZ and MC powers given in WXL, we have 

P„,/P„ = f (29) 
la 2—1 — sm B 

*' ^ screw \ ^/ 

and 

Pscrew/ Pq = 2a^ / /-, x • 2 /j " (30) 

2-(l-5)sin ^ 
where Pq = {B^hTM^ ^ BKM/Mq^ x 6.59 x lO^^ery • s'K 

By using equations (28) — (30) we obtain the curves of Pbz/Pq, Pmc/Pq ^-nd Pscrew /Pq varying 
with the parameters a* and n as shown in Figures 6 and 7, respectively. Defining the ratio of Pscrew to 
the total power Ptotai as r]p = Pscrew / Ptotai ■ wc have the curves of rjp varying with a* and n as shown in 
Figure 8. In the following calculation A = O.Stt is assumed. 

From Figures 6, 7 and 8 we obtain the following results: 

(1) For the given value of n, Pbz, Pmc and Pscrew all vary non-monotonically with the BH spin a*, 
attaining the maxima as a* approaching unity. 

(2) For the given value of a*, Pbz increases monotonically, Pmc decreases monotonically , while 
Pscrew varies non-monotonically with the power-law index n. 

(3) The average power Pscrew is always less than Pmc, and it is less than Pbz in most cases, except 
that the power-law index is small as shown in Figure 7a. 

(4) The ratio of power, rjp, increases monotonically with a*, while it varies non-monotonically with 
the power-law index n, attaining its maximum rjp ~ 0.084 at n ~ 5.25 with a* ~ 0.997. 

Therefore we conclude that the screw instabihty contributes only a small fraction of magnetic 
extraction of energy from a rotating BH. 



5. SUMMARY 

In this paper the screw instability of the magnetic field connecting a rotating BH with its surround- 
ing disk is discussed in the simplified model of CEBZMC. In our model the configuration of the poloidal 
component of the magnetic field is assumed, which remains constant on the horizon, and varies with the 
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radial coordinate on the disk according to a power law. By using an equivalent circuit for CZBZMC, 
we can calculate the poloidal current on the horizon and the disk, and then the toroidal component of 
the magnetic field on the disk is derived by using Ampere's law. Therefore, the criterion of the screw 
instability in our model is derived based on the Kruskal-Shafranov criterion. 

Although the discussion is carried out under some simplified assumptions, several interesting results 
for the screw instability have been obtained. 

(1) It is shown that the BH spin a* and the power- law index n are involved in the criterion, and 
the screw instability will occur if the two parameters are greater than some critical values. This result 
means that the screw instability is apt to take place for a disk containing a fast-spinning BH, and has 
the magnetic field concentrating in the inner region. 

(2) We find that the disk region for the screw instability can be also determined by a* and n by 
virtue of the criterion, and the screw instability is apt to take place outside some critical radius on the 
disk. 

(3) We prove that the state of CEBZMC always accompanies the screw instabihty as shown in the 
parameter space consisting of a* and n. 

(4) As an energy mechanism the effect of the screw instability on the magnetic extraction is dis- 
cussed. It turns out that the contribution of the screw instability is not significant in our model of 
CEBZMC. 

In this paper we fail to discuss the screw instability of the magnetic field in the BZ process, since 
we know very few about the remote astrophysical load. Considering stringent upper bound to the BZ 
power due to the screw instability in the BZ process (Li 2000a), we should give an overall evaluation 
on the effects of the screw instability on energy extraction from a rotating BH in the future work. 

Acknowledgements: This work is supported by the National Natural Science Foundation of 
China under Grant Numbers 10173004 and 10121503. 
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Fig. 1. — Poloidal magnetic field connecting a rotating BH with a remote astrophysical load and the 
surrounding disk. 
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Fig. 2. — Equivalent circuit for a model of CEBZMC. 
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Fig. 3. — Value of {2nzuD/L)F{a^;^,n) at ^ = 5 (a) varying with a* for n = 4 (b) varying with n for 
a* = 0.998. 
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Fig. 4. — Critical line for CEBZMC (thick solid line) and critical lines for screw instability with different 
values of factor A — O.Stt, O.Stt, 0.77r,7r in sohd, dashed, dotted and dot-dashed lines, respectively. 
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Fig. 5. — Contours of constant values of iscrew for the screw instability with A — O.Stt. 
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Fig. 6. — (a) Curves of Pb^/^o (solid line), Pmc/-Po (dashed line), and Pgcrew/Po (dotted line), and (b) 
curve of Pgcrew/ Pq (dotted line) varying with a* for n — h. 
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Fig. 7. — (a) Curves of Pb^/^o (solid line),PMc/-Po(dashed line), and Pgcrew/Po (dotted line), and 
curve of Pgcrew/ Pq (dotted line) varying with n for a* = 0.998. 
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Fig. 8. — Ratio of Pscrew to total power Ptotai (a) varying with a* for n = 5, (b) varying with n for 
a* = 0.998. 





